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Abstract 

We show that the effective diffusivity matrix D{V''^) for the heat operator dt — (A/2 — 
Vy"V) in a periodic potential V"" = ^22=0 Uk{x / Rk) obtained as a superposition of 
Holder-continuous periodic potentials Uk (of period T*^ := M'^/Z'^, d G N*, ^7^(0) = 0) 
decays exponentially fast with the number of scales when the scale-ratios Rk+i/Rk are 
bounded above and below. From this we deduce the anomalous slow behavior for a 
Brownian Motion in a potential obtained as a superposition of an infinite number of 
scales: dyt = doJt — W°°{yt)dt 
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1 Introduction 

Homogenization in the presence of a large number of spatial scales is both very impor- 
tant for applications and far from understood from a mathematical standpoint. In the 
asymptotic regime where the spatial scales separate, i.e. when the ratio between suc- 
cessive scales tends to infinity, multi-scale homogenization is now well understood. See 
for instance Bensoussan-Lions-Papanicolaou [111], or Avellaneda B, Allaire - Briane P], 



[31|, Jikov - Kozlov |^] or Avellaneda - Majda [|6[]. 

Nevertheless the case of multi-scale homogenization when spatial scales are not clearly 
separated, i.e. when the ratios between scales stay bounded, has been recognized as 
difficult and important. For instance, Avellaneda |^] (page 267) emphasizes that "the 
assumption of scale separation invoked in homogenization is not adequate for treating 
the most general problems of transport and diffusion in self-similar random media" . 
The potential use of multi-scale homogenization estimates for applications are numerous 
(see for instance |42| for applications to geology, or [^], 0] for applications to 



Differential Effective Medium Theories). The main application of this line of ideas is 
perhaps to proving super-diffusivity for turbulent diffusion: see for instance §, 0; H], 
II, m, 1^, m, ll, ||,|ll, H]; H, 0; H or I 



We are here interested in sub-diffusivity problems. Consider the Brownian motion in a 
periodic potential, i.e. the diffusion process 

dyt = doJt - VV{yt)dt (1.1) 

where V is periodic and smooth. It is a basic and simple fact of homogenization theory 
that yt behaves in large times like a Brownian motion slower than the Brownian Motion 
LOt driving the equation, i.e. y'^(t) = eyj/^a converges in law to a Brownian motion with 
diffusivity matrix Diy) < I^. 

We first treat here the case where ^ is a periodic n-scale potential with ratios (between 
successive scales) bounded uniformly on n.We introduce a new approach which enables 
us to show exponential decay of the effective diffusivity matrix when the number of 
spatial scales grows to infinity. 

From this exponential decay we will deduce the anomalous slow behavior of Brownian 
motions in potential V, when V is a superposition of an infinite number of scales. 
We have studied this question with a particular application in mind, i.e. to prove that one 
of the basic mechanisms of anomalous slow diffusion in complex media is the existence 
of a large number of spatial scales, without a clear separation between them. This 
phenomenon has been attested for very regular self-similar fractals (see Barlow and Bass 



llOl and Osada |^ for the Sierpinski carpet) (see also IQ). Our goal is to implement 
rigorously the idea that the key for the sub-diffusivity is a never-ending or perpetual 
homogenization phenomenon over an infinite number of scales, the point being that our 
model will not have any self similarity or local symmetry hypotheses. 
Our approach gives naturally much more detailed information in dimension one and this 
is the subject of []39| . 

This approach will be shown in forthcoming works to also give a proof of super-diffusive 
behavior for diffusion in some multi-scale divergence free fields (see jll[] for the simple 
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case of shear flow and [40| for a general situation). 

The second section contains the description of our model; the third one, the statement 
of our results and the fourth one the proofs. 



2 The multi-scale medium 

For U G L°°(T^) (we note := OT^), let mu be the probability measure on T'^ defined 
by 

mt/(dx) = e-2^(^)dx/ / e-2^(^)dx (2.1) 

The effective diffusivity D{U) is the symmetric positive definite matrix given by 

HD{U)1= inf / |/- V/(x)pm{/(dx) (2.2) 

for / in '§>'^~^ (the unit sphere of M'^). Our purpose in this work is to obtain quantitative 
estimates for the effective diffusivity matrix of multi-scale potentials V^" given by a sum 
of periodic functions with (geometrically) increasing periods: 

n 

Vo=T.^k{^) (2.3) 

fc=0 

In this formula we have two important ingredients: the potentials Uk and the scale 
parameters R^- We will now describe the hypothesis we make on these two items of our 
model. 

1. Hypotheses on the potentials Uk 
We will assume that 

Uk G C"(T'^) (2.4) 

Uk{0) = (2.5) 

Here C"(T'^) denotes the space of a-Holder continuous on the torus T'^, with < 
a < 1. We will also assume that the C"-norm of the Uk are uniformly bounded, 
i.e. 

Ka := supsup \Uk{x) - Uk{y)\/\x - y]" < oo (2.6) 
We will also need the notation 

i^o := supOsc([/fc) (2.7) 

fceN 

where the oscillation of Uk is given by Osc(C/) := supf/ — inf [/. 
We also assume that the effective diffusivity matrices of the Uk's are uniformly 
bounded. Let Xmin{D{Uk)) (respectively Xma.x{D{Uk))) be the smallest and largest 
eigenvalues of the effective diffusivity matrix D{Uk)- We will assume that 

Xrain- inf HD{Uk)l > (2.8) 
Amax:= sup HD{Uk)l<l (2.9) 
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2. Hypotheses on the scale parameters Rk 

Rk is a spatial scale parameter growing exponentially fast with A;, more precisely 
we will assume that Rq = = 1 and that the ratios between scales defined by 

rk = Rk/Rk-i G N* (2.10) 
for A; > 1, are integers uniformly bounded away from 1 and oo: we will denote by 

Pmin := inf rfc and pmax := sup (2.11) 

ken* ken* 

and assume that 

Pmin > 2 and Pmax < OO (2.12) 

As an example, we have illustrated in the figure the contour lines of V^^(x,y) = 
Efc=o^(^';^) ' ^ith /5 = 4 and U{x,y) = cos(x + 7rsin(y) + 1)^ sin(7rcos(x) - 2y + 
2) cos(7r sin(x) + y) 



3 Main results 

3.1 Quantitative estimates of the multi-scale effective dif- 
fusivity 

3.1.1 The central estimate 

Our first objective is to control the minimal and maximal eigenvalues of D{Vq'). More 
precisely writing 1^ the d x d identity matrix we will prove that 



Theorem 3.1. Under the hypotheses ^2. (\ ), i\2.H ) and p'^^^ > Ka there exists a constant 



C depending only on d, a, Ka, Kq such that for all n > I 

n n 

Id e-"^ n ^^in{D{Uk)) < D{Vo^) < h e"^ \{ \mUD{Uk)) (3.1) 

k=0 k=0 

where 

^ = Cp-t!^ (3.2) 
In particular e tends to 0, when p^i^ — > oo 

Remark 3.1. One can interpret this theorem as follows: /^(V^") is bounded from below 
and from above by the bounds given by reiterated homogenization under the assumption 
of complete separation between scales, i.e. pmin oo (product of minimal and maximal 
eigenvalues) times an error term e"*^ created by the interaction or overlap between the 
different scales. 

Remark 3.2. Originally the problem of estimating D{Vq) was called for in connection 
with applied sciences, and heuristic theories such as Differential Effective Medium theory 
have been developed for that purpose. This theory (DEM theory) models a two phase 
composite by incrementally adding inclusions of one phase to a background matrix of 
the other and then recomputing the new effective background material at each increment 



[|15[, [^], llTj]. It was first proposed by Bruggeman to compute the conductivity of a 



two-component composite structure formed by successive substitutions ( [|16| and ||l| ) and 
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Figure 1: A particular case 
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generalized by Norris |35f| to materials with more than two phases. 
More recently Avellaneda |Q] has given a rigorous interpretation of the equations obtained 
by DEM theories showing that they are homogeneous limit equations with two very 
important features: complete separation of scales and "dilution of phases". That is to 
say, each "phase" U^is present at an infinite number of scales in a homogeneous way. Yet 
two different phases never interact because they always appear at scales whose ratio is oo. 
Moreover the macroscopic influence of each phase is totally (but non uniformly) diluted 
in the infinite number of scales at which it appears. In our context, complete separation 
of scales would mean that R^+i/Rk grows sufficiently fast to oo, and "dilution of phases" 
would bmean that Fg" = J2k=o^ki^/^k) with UJ} — > as n ^ oo. The rigorous tool 



used by Avellaneda to obtain this interpretation is reiterated homogenization |12|. 



May be the most recent work on this topic is the article by Jikov-Kozlov who 
work under the assumption of "dilution of phases" and fast separation between scales, 
more precisely under the condition that J2T=i ^{^k/ Rk+i) < oo. Jikov-Kozlov use 
the classical toolbox of asymptotic expansion, plugging well chosen test functions in the 
cell problem. This method of asymptotic expansion is simply not at all available in our 
context. 



Theorem 3.1 will be proved by induction on the number of scales. The basic step in 
this induction is the following estimate ( |3.4| ) on the effective diffusivity for a two-scale 
periodic medium. 

Let U,T e C"(T°'). Let us define for R G N*, SrU G C"(T°') by SrU{x) = U{Rx). We 
will need to estimate D{SrU + T) the effective diffusivity for a two-scale medium when 
i? is a large integer. Let us define D(U,T), the symmetric definite positive matrix given 
by 

HD{U,T)1= inf f \l-Vf{x))D{U){l-Vf{x))mT{dx) for / G M'^ (3.3) 

/gC°=(Tf) Jxf 

Theorem 3.2. Let R e W and U,T e C"(T'^). // fl" > ||r||c, then there exists a 
constant C depending only on d,Osc{U), \\U\\a,ct such that 

e-'D{U, T) < D{SrU + T)< D{U, T)e' (3.4) 

with e = CR-"'/'^ 



Remark 3.3. Theorem 3.2 implies obviously that 



D{U,T)= \iui D{SrU + T) (3.5) 

R—*oo 

SO that D{U, T) should be interpreted as the effective diffusivity of the two scale medium 
for a complete separation of scales. Naturally D{U, T) is also computable from an explicit 
cell problem (see ( [4.12D ). 



Remark 3.4. The estimate given in theorem |3.2| is stronger than needed for theorem 3.1 



It gives a control of D{SrU + T) in terms of D{U, T) and not only of the minimal and 



maximal eigenvalues of D{U) and D{T). In fact we will only use its corollary 3.1 given 



below, which is deduced using the variational formulation (3.3) 



Corollary 3.1. Let G N* and U,T e C"'(T^). // i?" > ||r||c, then there exists a 
constant C depending only on d,Osc{U), \\U\\a,Oi such that 

\n.in{D{U))D{T)e-' < D{SrU + T) < X^,^{D{U))D{T)e' (3.6) 
with e = CR-""'^ 
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Remark 3.5. We mentioned that theorem 3.1 is proved by induction. This induction 



differs from the one used in reiterated homogenization or DEM theories by the fact we 
homogenize on the larger scales first and add at each step a smaller scale. 

Let us introduce the following upper and lower exponential rates 



Definition 3.1. 



A+ = lim sup - In A^ax(l?(V^o")) (3-7) 

n— »oo ri 



A" = lim inf - In Amm(-D(K")) (3- 

n^oo n 



Theorem |3.1| implies the exponential decay of D{Vf^ i.e. 
Corollary 3.2. Under the hypotheses ( \2.6[ ), { 2.1O0 and p'^^^ > Ka one has (with e given 



by forn>l 

he-^^K^l < D{V,n < he-^KLi (3.9) 

and 

A+<lnAr„ax + e (3.10) 



A">lnA„,in-e (3.11) 

In particular if Xmax < 1 then there exists a constant po = (l + Cd^Ko,Ka,a/ {— ''^max))'^ 
such that, for p^am > Po 

A+ < (3.12) 

Thus one obtains the exponential decay of D(V^) only for a minimal separation 
between scales, i.e. Pmin greater than a constant pQ characterized by the medium. It 
is natural to wonder whether this condition is necessary and what happens below this 
constant pQ. We will give a partial answer to that question, in the simple case when the 
medium V is self-similar. We will see that it is possible to find models such that for a 
certain value C of the separation parameter pmin = C, D{Vf^) decays exponentially and 
for Pmin = C + 1, D(y^) stays bounded away from zero. This will be done using a link 
with large deviation theory. 

3.1.2 The self similar case 

Definition 3.2. The medium V is called self similar if and only if Vn, Un = U and 
i?„ = with pen, p>2 

Definition 3.3. For U € C"(T'^) and p G N/{0, 1} we denote by pp{U) the pressure 
associated to the shift Sp{x) = px on T'^, i.e. 

PpiU) = sup ( / [/(x)d^(x) + hpip)) (3.13) 

where hp is the Komogorov-Sinai entropy related to the shift Sp. We denote 

PpiJJ) = pp{U) - pp{0) = Pp{U) -dlnR (3.14) 
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We refer to [^] and [41| for a reminder on the pressure, let us observe that Vp{0) 
differs from the standard definition of the topological pressure by the constant din p so 
that VpiO) = 0. 

We will relate in the self-similar case the exponential rates A"^ and to pressures for 
the shift Sp and to large deviation at level 3 for i.i.d. random variables. 
In the self similar case we will write \~{—U) the exponential rates associated to D{—Vf^). 
We write 



Z{U) 



Vp{2U)+Vp{-2U) 



(3.15) 



Theorem 3.3. If the medium V is self similar then 
1. Ifd = l 

X+(U) = x-{u) = Z{U) 



(3.16) 



2. ifd = 2 then 



X+{U) + X-{-U) = Z{U) 



(3.17) 



Moreover if there exists an isometry A of such that U{Ax) = —U{x) and a 
reflection B such that U{Bx) = U{x) then A~(-C/) = \'{U) = A+(C/) so that 



X+{U) = X-{U) = Z{U)/2 



(3.18) 



3. For any d 



Z{U) < X-{U) < A+([7) < 



(3.19) 



Remark 3.6. The statement ( 3.16 ) is obtained from the explicit formula for D(VJ)"') in 
d=l, see [HI. 



Remark 3.7. It is obviously important to know when Z{U) is strictly negative to be able 
to use this theorem. A well known and useful criterion can be stated as Z{U) < if and 
only if U does not belong to the closure of the vector space spanned by cocycles, which 
can be shown to be equivalent to say that 



Z{U) < 4^ lim sup -\ V {U{p^x) - [ U{x)d. 

n^oo n II Jjd 



> 



(3.20) 



We refer to |39| for the proof of the last statement. 



Example 3.1. Let U{x) = sin(x) — sin(81x) in dimension one. In fact (I3.16D and (^^ 
shows that A"^([/) < as soon as p > 82. For p < 81 the situation is a bit surprising: 
A+(C/) < for p / 3,9,81. For these exceptional values A+(C/) = and in fact D{V^) 
remains lower bounded by a strictly positive constant. 

This example shows that for a given potential U , even though the multi-scale effective 
diffusivity Dly^) decays exponentially for p large enough, one can find isolated values 
of the scale parameter for which Dly^) remains bounded from below. 

I are only used to prove 
D{Vq)) (see proposition 



Remark 3.8. The symmetry hypotheses given in theorem 
that for all n, D{V^) = D{-V^) and X^^^{D{V^)) = X^m 

m 
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3.2 Sub-diffusive behavior from homogenization on infinitely 
many scales 

Here we consider the diffusion process given by the Brownian Motion in the potential 

oo 

^ = ^o°° = 2^f^fc(^M.) (3.21) 

fe=0 

We assume in this section that the hypotheses (U) , (U), (U), (HJ), and 

( 2.12| ) hold. To start with we will assume that 

a = 1 and that the potentials C/^ are uniformly C^. (3.22) 

In particular V is well defined and belongs to C"'^(IR'^) and ||Vy||oo < oo. 

The diffusion process associated to the potential V is well defined by the Stochastic 

Differential Equation 

dyt = diot - VV{yt)dt (3.23) 

We will show that the multi-scale structure of V can lead to an anomalous slow behavior 
for the process yt- To describe this sub-diffusive phenomenon we choose to compute the 
mean exit time from large balls, i.e. Let 

r(r) = inf{t > : \yt\>r} (3.24) 

We would like to show that E2.[T(r)] grows faster than quadratic in r when r — > oo 
uniformly in x. We cannot obtain such pointwise results in dimension d > 1 (see 3.2.1 
for a discussion, the case d = 1 is treated in |^9|). But we will start with averaged results 
on those mean exit times. 

The fact that the homogenization results of subsection p.l| can be of some help to estimate 
the mean exit times is shown by the following lemma 

Lemma 3.1. For U E C°°{Tj^), (R > 0) writing the exit times associated to the 
diffusion generated by Lu = A/2 — VC/V one has 
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(3.25) 



.2 



Let mv,r be the probability measure on the ball B{0,r) given by 

mvAdx) = e-2^(^) dx/{ [ e'^^^^') dx) (3.26) 

' JB{0,r) 

We will consider the mean exit time for the process started with initial distribution mv^r, 
i.e. 

^my,[r{r)]= [ E^[T{r)]mvAdx) (3.27) 

JB(0,r) 
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Theorem 3.4. Under the hypothesis Amax < 1 there exists C2 depending on d, Amax; Kq, Ka,oi 
such that if Pmin > C2 then 

lim inf "^■'•^ ^ > 2 (3.28) 

More precisely there exists C3 > 0, C4 > 0, C5 > such that for r > C3, 

^my,M^)] =r2+'^(^) (3.29) 

with 

< Q < }^^S(i _ _ < iy{r) < 5^(1 + ^1^) + J-C^ (3.30) 

ill Pmax ill Pmin 

Inr mpmin inpmin mr 

Where C3 and C5 depend on (d, Kq, Ka, a) and C4 on (Amax, Pmax)- 



The proof of this result rehes heavily on theorem 3T. The idea being that ^mv,- [''"(^)] 
is close, when r is large to r^/ Amax (-^^(^o")) where n is roughly sup{m € N : Rm l£ r}. So 
that the exponential decay of D{V(^) gives the super-quadratic behavior of ^mvr [''"(^)] ; 
i.e. sub-diffusivity. 



Remark 3.9. The differentiability hypothesis ( 3.22| ) though convenient in order to define 



the process yt as a solution of the SDE |3.23 is in fact useless. The theorem is also 



meaningful and true with < a < 1. See section 4.2 for an explanation 



3.2.1 Pointwise estimates on the anomaly 



Theorem |3.4| gives the anomalous behavior of the exit times with respect to the invariant 
measure of the diffusion and it is desirable to seek for pointwise estimates of this anomaly. 
The additional difficulty is to obtain quantitative estimates on the stability of divergence 
form elliptic operators under a perturbation of their principal parts (see conjecture |3.1| ). 
By stability we mean here the validity of the following condition |3.l|. 



For U € C^{B{z,r)). Write E^, the expectation associated to the diffusions generated 
by Ljj = — VC/V. V is said to satisfy the stability condition if and only if: 



Condition 3.1. There exists fi > such that for all n £ N, all z G M'^, and all r > 0, 
le-MOscs(,,,)(V-i) -j^f E^°[T{B{z,r))]<EY[TiB{z,r))] (3.31) 



^]"[r(S(z,r))] </ie^O'^=s(.,.)(v;°Ji) g^p [T{B{z,r))] (3.32) 

xeB{z,r) 

Where OscB(z,r){U) stands for sup^^^.r) ^ ~ '^^^B{z,r) U . Under the stability condition 



3.1, we can obtain sharp pointwize estimates on the mean exit times. 



Theorem 3.5. If V satisfies the stability condition 3J_, then there exist a constant Cq 

depending on {d, KQ,Ka, a, fj,, Amax) such that for Pmin > Cq, one has for all x € M"^ 

InE J T(B(x,r))] 

lim inf \^ ^ > 2 (3.33) 

r— »oo mr 

More precisely there exists a function a{r) such that for r > one has 

^g^2+a{r){i-7) < [t{B{x, r))] < C79r2+'^('^)(i+^) (3.34) 
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with 

+ <f7(r) < + (3.35) 

^max ill Pmin '-^ Pvain '-^ Pmin 

and 7 = Cs-fCo/Clnpmm) < 0.5. Where the constants C^,C/i,C-j,C% and Cg depend on 
(d, Ko,^Q,a,/", Amax) and on d. 

Here T{B{x,r)) denotes the exit time from the ball B{x,r). 
Remark 3.10. In fact cr(r) can be described rather precisely. Let, 

a(r, n) = - In \^^^D{V^) j \nr (3.36) 

Define nef{r,Ci,C2) = sup{n > : e^^+^^^^^oi^^ < ^^^2^ n^^ien there exists Ci,C2 



depending only on d such that (T(r) in theorem 3^ is a{r, ne/(r, Ci, C2)). 



Using the precise information of theorem 3^ we can estimate the tails of probability 
transitions for the process yt (or the tail of the heat kernel for the operator Ly). We get 
non-Gaussian upper bound similar to the (more precise) ones proved for fractal diffusions 
(see and [||]) 



Theorem 3.6. IfY satisfies the stability condition 3_A, then for Pmin 

> Ce, r > 

Cwr < t < Cnr2+'^('^)(i-3^) (3.37) 

one has 



InF^U -x\>r]< -CisjC-y^'^''^ (3-38) 



with (Cn < O.SlnpminJ 



< I" (1 - -^^) < u{y) < |^!^^(1 - -^^) (3.39) 
Where Cq,Ci3,Ci4 and C15 depend on (d, Kq, Ka, 

(d, Kq, Cy, /i, Ajnax) Pmim Pmax) 



Remark 3.11. The non Gaussian structure of ( 3.38 ) is similar to the one obtained for 
diffusion processes in fractals. Indeed 

-C^y(;^r = -C(^^)^ (3.40) 

with ~ 2 + I/. 

Next, it has been shown in [g9[ that for U G L°°(T^), lnp^(x,y,t) is roughly — *(y — 
x)D~^{U){y — x)/t for t > — y| (homogenized behavior). Where p^ is the heat kernel 
associated to Lu. Next, writing nef{t/h) = sup„{ii„ < t/h} the number of scales that 
one can consider as homogenized in the estimation of the heat kernel tail one obtains 
from a heuristic computation (which can be made rigorous in dimension one, see |39]) 
that for Cioh<t< Cuh^+'', 

lnP(y, >h)< -C^^^^ r. -C-{-r^ r. (3.41) 
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with dyj ^ 2 — The equation 3.41 suggests that the origin of the anomalous shape 
of the heat kernel for the reflected Brownian Motion on the Sierpinski carpet can be ex- 
plained by a perpetual homogenization phenomenon and the formula linking the number 
of effective scales and the ratio t/h. 

The condition Cioh < t can be translated into "homogenization has started on at least 
the first scale" (ng/ > 1) and the second one t < Cuh?^^ into "the heat kernel associated 
to Ly is far from its diagonal regime" (one can have h? /t << 1 before reaching that 
regime, this is explained by the slow down of the diffusion). 



The weak point of theorems 3.5 and 3.6 is naturally that checking condition 3.1 
seems difficult. But we believe that in fact this condition is always true (we refer to the 
chapter 13 of (3^]), since this condition is a consequence of the following conjecture (see 
[l39| , proposition 2.3). 



Conjecture 3.1. There exists a constant Cd depending only on the dimension of the 
space such that for A G C°°(5(0,l)) such that A > on 5(0,1) and 4>,ip £ (7^(5(0,1)) 
null on dB(0,l) and both sub harmonic with respect to the operator — V(AV), one has 



(3.42) 



X{x)\V(j){x).Vij{x)\dx <Cd X{x)V(j){x).Vi;{x)dx 

B(0,1) Jb{0,1) 



It is simple to see [J 
proving conjecture |3T 
estimate for the heat kernel in theorem 3.(: 



that this conjecture is true in dimension one with = 3. So 



would give the pointwise estimates of theorem 3.5 and the tail 



Remark 3.12. Here we have assumed that the are uniformly Ci but let us observe 
that since theorems p7 



and 3^ are robust in their dependence on a and (one can 
choose a < 1). One can built a process, with the assumption that the Uk are o-Holder 
continuous, whose mean exit times and heat kernel tail verify the estimates given in 



theorems 3.5 and 3.6. 



4 Proofs 

4.1 Multi-scale homogenization with bounded ratios 

4.1.1 Global estimates of the multi-scale effective diffusivity: theorem 



3.1 



The proof of theorem 3.1 will follow from the Corollary 3.1 by a simple induction. Let 



n G A^/{0, 1}, p G N, 1 < p < n and assume that 

n n 

j,e-(n-p)e(p_) -Q A^i,(Z)(C/,)) < D{V^) < I^eS^-PHP^-^'^) J] Xm.AD{Uk)) (4.1) 

k=p k=p 

One pass from the quantitative control on D{Vj!^) to a control on D{V^_i) by choosing 
U{x) = Up-i{x), T{x) = Vp{Rnx) and R = Rn/Rp^i in theorem ^ and observing that 
|[r||a/i?° < (2" — l)^^Ka/p'^J^^^^■ This provcs the induction and henceforth the theorem. 



4.1.2 Quantitative multi-scale-homogenization: Upper bound in the 
theorem |3.2| 



4.1.2.1 We will use the notation introduced in theorem 3.2. By the variational 
formula (^), D{U) is continuous with respect to U in L°°-norm, it is sufficient to prove 
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theorem |3.2| assmning that U and T are smooth. 

First let us prove that when homogenization takes place on two scales separated by a 
ratio R, the influence of a translation of the first one with respect to the second one on 
the global effective diffusivity can easily be controlled, i.e. for y G T{^, writing Qy the 
translation operator T{x) QyT{x) = T{x + y) 

Lemma 4.1. 

e-^^^D{SRU + T)< D{SnQyU + T)< e^^^DiSnU + T) (4.2) 

Proof. The proof follows by observing that SrU + QyT = @[Ry]/R{SiiU + T) + QyT — 
where [Ry] stands for the vector with the integral parts of {yR)i as coordinates. 
Thus by the variational definition of the effective diffusivity 

D{SrU + OyT) < e^ll®«^-Q[««i/«^ll-Z)(e[Rj^]/^(5/j?7 + T)) (4.3) 

and the equation ( [4. 2D follows by observing that the effective diffusivity is invariant under 
a translation of the medium: D{@[RyyR{SjiU + T)) = D{SrU + T). □ 

Next we will obtain a quantitative control on Jy^jd D{SrU + QyT)dy 

Lemma 4.2. For R > ||r||Q 

[ D{SrU + @yT)dy <e^^^-^{l + Cde^''^'''^^\\\T\\a/R'')^2)D{U,T) (4.4) 



Let us observe that the combination of lemma 4.2 with 4.1 gives the upper bound 



in theorem 3.2. 

Write xY the solution of the cell problem associated to U. We remind that for / € 
Lu = 1/2A - Vf/V, Luxi = -IVU, xYiO) = and 



HD{U)1= / \l-VxiVmu{dx)= / \l-Vxi).lmu{dx) (4.5) 

Write 'x^'y^^''^ the T*^ periodic solution of the following cell problem (which corresponds 
to a complete homogenization on the smaller scale): for / € S''"^ 

V(e-2^i?(C/)(/-Vxf(^^''^))=0 (4.6) 

Write for y E T'^, x x{x, y) the solution of the cell problem associated to SRQyU + T. 
Let / G S*^"^, by the formula associating the effective diffusivity and the solution of the 
cell problem and using that / — VxXi{x,y) is harmonic with respect to Lsj^u+OyT, one 
obtains 

[ HDiSRQyU + T)ldy= [ {1-Vxxi{x,y)).ldxdy (4.7) 
Writing the decomposition 

l = {Id- Vx^'iRx + y)){l - Vxf ^""^'^(x)) + Vx^(fe + y)(Vxf ^""^'^(x) - l)mu(^R.:,+y)+Ti.){dx) dy 

(4.8) 

we get that 

/ HDiSRQyU + T)ldy = h-l2 (4.9) 
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With 



{I - ^xXi{x,y)){Id - Vx^ {Rx + y)){l - '^xf"'^^''^ {x))mu(R.+y)+T{.){dx) dy 

(4.10) 



and 



h= (l -'^xXi{x,y)Wx^{Rx + y){Vxf''-^^''^{x) -l)mu,^R^+y^j^T{.)idx)dy 

(4.11) 

It is easy to see that x^^^^''^ is a minimizer in the variational formula (3.3) associated 
to D{U,T), which is the effective diffusivity corresponding to two-scale homogenization 
on U, T with complete separation between the scales, that is to say: 

D{U,T)= [ *(I,-Vx^(^)'^(x))Z?(C/)(/rf-Vx^(^)'^(x))mT(dx) (4.12) 

A simple use of Cauchy-Schwarz inequality gives an upper bound on Ii, 



(x',j/)G{T<*)2 
(a;,j/)g{T'*)2 



\l - VxXi{x,y)?m^^''^+y^+^^^\dx)dy 



\{h - Vx'/iRx + y)){l - Vxf ^'')'^(x))|2m^(^-+S')+^(-)(dx), 



(4.13) 



Integrating first in y in the second term and, using the formulas linking effective diffu- 
sivities and solutions of the cell problem, we obtain 



h < 



yeT^ 



lD{SRQyU + T)ldy) X riD{U,T)l) e—' 



(4.14) 



We now estimate I2. The following lemma together with (4.9) and ( 4.14| ) gives lemma 



4.2 



Lemma 4.3. 



\h\< 



HD{SRQyU + T)ldy) X (HD{U,T)1 

The proof of this lemma relies heavily on the following elliptic type estimate 
Lemma 4.4. 

Ilxflloo <Crfexp((3d + 2)Osc(i7))|Z| 



(4.15) 



(4.16) 



This lemma is a consequence of theorem 5.4, chapter 5 of |43] on elliptic equations 
with discontinuous coefficients (see also [0]), we give the proof of lemma 4.4 for the sake 
of completeness in paragraph [4.1.2.2 . 

We will now prove lemma |4^ . First we will estimate the distance between xi{x,y) and 
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Xi{x + y / R, 0) for y G [0, l]'^ in norm. 

By the orthogonality property of the solution of the cell problem for y G [0, l]*^ 

f y ^-2(U{Rx+y)+T(x)) 

-2{U{Rx+y)+T(x)) 



(4.17) 



<HD{SbU + T)/e^^ - HD{SRQyU + T)l 



Thus by lemma for y G [0, 1] , 

WxXi{x+ ^,0) - Va;X«(3;,y)pnT,[;(^,^+j^)+2^(.)(dx)(iy 
x&Td ^ (4.18) 

< HDiSnQyU + T)l{e^^ - 1) 
Let us introduce 



h= [ (Z-V,xKx + ^,0))Vx^(i?x + y)(Vxf^^^'^(x)-0 



(x,j/)eT<*x[0,l] 

^)+T{x+^)) 



^-2iUiRx+y)+T{x+^)) 



(4.19) 



one has 



Or - - 

HDiSROyU + T)l dy) ' X (HD{U, T)l) ' 
yeTd ^ ^ ^ (4.20) 

gg4Q^+Osc(l/)(g8&^ 

This can be seen using ( |4.18D , ( |4.11| ) and by Cauchy Schwarz inequality (the computation 
is similar to the one in (l4.13| )); Voigt-Reiss' inequality {D{U) > g-20sc{c/)'j noticing 
that for y G [0, l]'^,\T{x + y/R) - T{x)\ < \\T\\^/Ra. 
We now want to estimate I3. Noting that 

V,(e-2(^(«^+^)+^(^+*))(/ - V.x/(x + |, 0))) = 

Vx^(iix + y) = Vj,x^(i?x + y) 
and integrating by parts in y, one obtains 

/3=E/ (e--(-*)(^-V.xKx + ^,0)) 



i=l 

e 



2T{.+^)(; _ v.xKx + |, 0))) .ea^(i?x + y^) 

D(U) T (,~2U(Rx+t) 

Where we have used the notation 

a^([0,l]'^) ={xG [0,lf : X, = 0} 
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Let us introduce 

d 



eTrf,j/»e9'([o,i]d) 



(4.21) 



It is easy to obtain 
I/4-I3I < de3 0^^(^)eT(eT - l)('*/D(5^jC/ + r)/ 



ID{U,T)1]' (4.22 



'IIy^II ^* 

Ha. IIoo 



We will now put into evidence the fact that although xK^j 0) '^ot periodic on R~^T'^ the 
distance (with respect to the natural norm) between the solution of the cell problem 
Xi(X)O) s-iid its translation Xl{x + ^k/R-,^) along the axis of the torus R'^T"^, is small. 
This is due to the presence of a fast period R~^T'^ in the decomposition V = SrU + T. 
Using the standard property of the solution of the cell problem one obtains 



'^Xi{x+ %,0) - Vxi{x,{))\^msj^u+T{dx) 

|/ - Vxi{x, 0) + Vxi{x + 5,0)- Vxi{x, Q)\^ms^u+T{dx) - HD(SrU + T)l 



< e fl° 



1/ - Vxz(x + ^, O))pm0^(5^^+^)(dx) - */Z)(5fiC/ + T)l 



(4.23) 



which leads to 



1 |VxK^ + - Vx/(x,0)|2m^«^+^((ix) < HD{SRU + T)l{e'''^ - 1) (4.24) 

Combining this inequality with the definition ( 4.21| ) of I4, Cauchy-Schwarz inequality 
proves that 



l^4|<E 



{{-yxXi{x + 



y' + ej 
R 



,0) 



X^(i2x + y*)(Vx; 



D(C/),T, 



0)^ 



g-2l/(/ir+j/0-2r(x+2j) 

J^,e-mz)dzJ^,e-'^nz)dz' 



dx ] dy^ 



Combining this with ( [4.24 ) one obtains 



h\ < {e^ -l)-2riD{SRU + T)l]U\\x 



„3 0sc(l/) ( t 



ID{U,T)1] e— 



(4.25) 



(4.26) 



Using lemma [4. IE to estimate ||x ||cx3 in (|]26|) and combining ( |422D and (|]2^ one 
obtains ( 4.15| ) and lemma O, which proves the upper bound of theorem |3.2| . 



16 



4.1.2.2 The purpose of this paragraph is to prove the estimate ( 4.16 ). First we will 
remind a theorem concerning elhptic equation with discontinuous coefficient from G. 
Stampacchia. Its proof in a more general form can be found in [^ ], chapter 5, theorem 
5.4 (see also Q). 

Let us consider the operator (in the weak sense)-L = V(AV) defined on some open set 

n C M'^ (for d > 3) with smooth boundary A is a d x d matrix with bounded 

coefficients in L°°(17) such that for all C e K'^, A|CP < ^^A^ and for all ij; \Aij\ < M, 

for some positive constant < A, Af < oo. 

Let p > d > 3. For 1 < i < d let fi G LP{n) 

if X £ -f^/oc(^) is ^ local (weak) solution of the equation 



V{AVx) = -^da, 



(4.27) 



i=l 



then X is in L°^{Q) and ii xq G Q and R > 

Theorem 4.1. The solution of ( 4-2'^) verify the following inequality (in the essential 
supremum sense with Vl{xq, R) = 0,0 B(xq, R) ) 



max |x[ < K 



n(xo,|) 



R"^ 



\X\\ 



n{xo,R) 



i=l 



1 — s 

R^ p 

i\\LP{n{xo,R)) — — 



(4.28) 



with K = Cdif)- 

The explicit dependence of the constants in M and A have been obtained by following 
the proof of G. Stampacchia |44|. We will now prove (P^ ) for d > 3 (For d = 1, this 
estimate is trivial, for d = 2, it is sufficient consider U{xi,X2) as a function on Tf to 
obtain the result), xi satisfies 

V{exp{-2U)Vxi) = /.Vexp(-2C/) 

then by theorem 4J. for xq G [0, l]'^ 

max Ixil < C'dexp(3 0sc(C/)(i)[( / [x^P) ^ + I^I exp(2 Osc(f/))] 
Now by periodicity 



Ixil'^dx < / \xi?dx 

B{xo,l) 



and by Poincare inequality (we assume Jja xi{x)dx = 0) 

/ \xi?dx<Cd [ \Vxi?dx 



thus 



And since 



one has 



/ \xi\^dx<Cdexp{2 0sciU)) [ \Vxi\^mu{dx) 

Jb{xo,1) JT<i 



fd 



\l-Vxi?mu{dx) = f 



fd 



\'^Xl?'mu{dx) 



j 



\Vxi\^muidx) < f 



and the bound on ||x/||oo is proven. 
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4.1.3 Quantitative multi-scale- homogenization: Lower bound in the 
theorem |3.2| 

4.1.3.1 As for the upper bound it is sufficient to prove theorem |3.2| assuming that 



U and T are smooth and we wiU use the notation introduced in the paragraph 4.1.2.1 . 
We wiU prove below that 

Lemma 4.5. If R > \\T\\a then for ^ G S'^-^ 

/ '^D{SrU + QyTY^^dy < (1 + Q,o.c(c;),||c/||.,a,||T|U^-"/')*e^(C^, T)-^i (4-29) 



This lemma with lemma [4.1| gives the lower bound in theorem We now prove 



lemma 4.5. Let us introduce 



Jjaex.p{-2{Sii@yU + T){x))dx 



(4.30) 

^"f^) - - J^^&jfk-J'^ - (4.31) 

and 

P^(^).^(x) = J, - j-^-^^D{U){I, - Vx^(^)'^(x).)Z?(C/,r)-^ (4.32) 

We remind that P{x, y) minimize the well known variational formula associated to 
D{SR@yU + r)~^, that is why it will play for the lower bound in the theorem |3.2| the 
role played by the gradient of the solution of the solution of the cell problem Vx{x,y). 
for the upper bound. More precisely, for ^ € one obtains as in the proof of the 

upper bound (by decomposing ^ here) 

/ 'iDiSaQyU + T)-'idy 

= / {I e-2(^«0^^+^)(^)dz)e2(5«e-^+^)(-)*e*(Id-P(x,y))eda;dy 

(4.33) 

with 

I,= [ e~2^(^) dz [ e-2^(^) dz [ ^2iSneyU+T)i.)t^t^j^ _ p^^^ 

JT'i Jj'i J(x,y)G(T'*)2 (4.34) 

{Id - P^'iRx + y)){Ia - P^(^)'^(x))C dxdy 



and 



/2=/ e-^^^'Uz [ e'^^^'Uz [ e2(^«®«^+^)(")*C*(/rf-P(x,y)) 

Jt'^ Jfl J (x,y)£{T'i)'2 (4.35) 

P^{Rx + y){Id- P^^^^'^{x))^dxdy 
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As for the upper bound, using Cauchy Schwarz inequality for the integration in x and 
y, and using 

^(D{U,T)~'^^= f ^2{SReyU+T)(x) 
((/, - P^'iRx + y)){h - 

one obtains that 



(4.36) 



|/i I < e 



'iD{SRQyU + T)-^i dy) ' ('^D{U, T)-'0 ' 



(4.37) 



Thus I2 will be an error term and it will be proven below that 
Lemma 4.6. 



\h\ < 



d,Osc{U),\\U\\a,at 



(4.38) 



1)^ 



Let us observe that combining the estimate ( 4.38[ ) of lemma 4.6 with ( [4.37| ) and ( f4.33 ) 
proves lemma |4.5|. 



We will now prove the lemma 4.6. As it has been done in the proof of the upper bound 
it is easy to show that, with 



y 



^2{U{Rx+y)+T(x+l)) 



(z,i/)gT<*x[0,1]<' 



(4.39) 



'e{h - P{x + -|, 0))P^(i?x + y){h - P^(^)'^(x))^ dxdy 
one has 

1/3 - /2I < 6e^'<^h^ (e^^ - 1)^ ('CDiSRU + T)A ' Nd{U, T)^ 



(4.40) 



It will be proven in 4.1.3.2 that 

Lemma 4.7. There exists dxdx d tensors such that Hf-^ = -H^^^ £ C°°(T^), 

d 

PjL = '^jH^'j^ and ||-H'(^m||oo < Q,Osc(;7),||(7||a,a (4-41) 



Combining ( |4.41|) with the explicit formula (4.30) for P one obtains 



/Tdexp(-2(5RC/ + r)(z))dz y(^,y)eTrdx[o,i]^.j5/ ^ ^^"^ ^ (4.42) 
dkHY^^^^iRx + y)({h - Vx.(x + ^))D{SrU + Py^i) . dxdy 

Thus, using the same notation as in the equation ( 4.21| ) and integrating by parts in y, 
one obtains 



E 



fja exp(-2(5ijC/ + T){z))dz . J{x,y'')€T'ixd''ilo,i]'') 



(*C*(/,-P^(^)'^(x))). 



i/{^,/i?x + /) ((Vx.(x + ^) - Vx.(x + ^^^))DiSRU + T)-'C) dx dy 



(4.43) 
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Which, using Cauchy Schwarz inequahty, leads to 

d 



ijk . _-, iy'=ea'-'([o,i]'^) ^ JxeTt* 



„2(J/(fc+3/'=)+T(x+iL))^^ 



fc=i^y'=ea'--([o,i]'^) 

((V7 f V7 / , + (4.44) 
((Vx.(x + — ) - Vx.(2; + 5 — )) 



D{SrU + T)-^C) e-2(^(«-+?'')+^(-+ ) ' dy^ 



Using bounds ( [4.41| ) and the equation ( |4.24D to control the natural distance between 
the solution of the cell problem x.(^ + U^/R) and its translation by e^/i? one obtains 



Combining (4.45) and ( 4.40| ) one obtains ( |4.38| ), which proves lemma 4.38 



(4.45) 



4.1.3.2 In this paragraph, we will prove lemma 4.7. Since for each m € {1, . . . , d}, 
are divergence free vectors with mean with respect to Lebesgue measure, by 
the proposition 4.1 of |28| there exists a skew-symmetric T^-periodic smooth matrices 
^^^P • • • > {ii'^.m = -Hf,J such that for all m 



U 

ijm 



Moreover writing 

pU = V T)'' p2i7T{k.x) 

the Fourier series expansion of P^, one has (see the proposition 4.1 of ||2 

^k u . „k 



^ _}_ Sr^ Pnmkj Pjm^n oiTrjk.x) 



Let us observe that 



(4.46) 



(4.47) 



(4.48) 



H'^. = - P" 



(4.49) 



where Bnm are the smooth T'^-periodic solutions of ABnm = djP^^. By theorem iA (the- 
orem 5.4, chapter 5 of ^^), if Bnm is chosen so that /^^ Bnm{x)dx = then ||-Bnm||oo < 

-^nmlloo- Now usi ng the orem 1.1 of [^] it is easy to obtain that ||VxPl|co ^ C'd,Osc(i/),||(7||a,a' 
combining this with ( 4.31 ), one obtains 



|-B;^mlloo < C'd,Osc((7),||t/||c«,o 



(4.50) 



Which leads to (|4.41j ) by the equation (|4.4g| ). 
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4.1.4 Explicit formulas of effective diffusivities from level-3 large devi- 



ations. Proof of theorem 3.3 



The equation ( |3.19| ) follows from the Voigt-Reiss inequality: for U G L°°{T'^) 

D{U)>h{[ e2^(^) / e-2t/W)-i (4.51) 

and the fact that, if ?7 G C"(T'^), then by the Varadhan's lemma and Ievel-3 large 
deviation associated to the shift Sp, 

\im-ln([ e^"^-oUi4^)dx) = PJU) (4.52) 

n^oo n \ Jfd ) 



We refer to |39| for a more detailed proof of this statement. 

In higher dimensions, when the medium is self-similar one can use the criterion ( p. 20 ) 
associated with the equation ( |3.19| ) to characterize ratios for which D(VJ)") does not 
converge to with an exponential rate. The equation ( ^.17 ), i.e. the extension of the 
result (3.16) to dimension 2 is done by observing that 

Proposition 4.1. For d = 2 one has 

Amax(l?(f/))A^in(l)(-C/)) = Knn{DiU)) X^^,{D{-U)) 

^ 1 (4.53) 

Jrpd exp{2U)dx Jrpd exp{—2U)dx 

From which one deduces that if D{U) = D{—U) then 

\unn{D{U))X^^^{D{U)) = { f eM'^U)dx [ eM-^U)dxy^ (4.54) 

Let us observe that the assumption D{U) = D{—U) is satisfied if, for instance —Un{x) = 
Un{—x) or —Un{x) = Un{Ax) where A is an isometry of W^. And the existence of a 
reflection B such that U{Bx) = U{x) ensures that Amm(-D(C/)) = \ma.^{D{U)). Thus 
these symmetry hypotheses combined with ( 3.17| ) ensure the validity of ( |3.18| ). 
It would be interesting to extend the equation ( |3.17 ) of theorem |3.3| to more general 
cases and higher dimensions. Indeed the proposition 4J is deduced from the following 
proposition |4.2| that put into evidence a strong geometrical link between cohomology and 
homogenization . 

Write Tsoi = {p e {C'=^{T'})Y\di\{p) = and Jj,dpdx = o| and Q{U) the positive, 
definite, symmetric matrix associated to the following variational problem. For / € 

[rpd \l - exp{2U)dx 
HQ{U)1 = inf — ^ (4.55) 

p<^^soi jrpd exp{2U)dx 

Write in the increasing order X{D{U))i and decreasing order X{Q{U))i the eigenvalues 
oiD{U) and Q{U). 

Proposition 4.2. For a// i G {1, . . . , d} 

X{D{U))A{Qm^ = 7 T^JTKT^r ( 9rrw (^"^^^ 

Jrpd exp[2U )dx Jj,dexp{—2u )dx 
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Now we will introduce a geometric interpretation of homogenization that will allow 
us to prove proposition 4.1 and equation (3.17) of theorem 3.3. Let U G C°°{T'^). It is 



easy to obtain the following orthogonal decomposition 

H = {L\mu)f = Hpot e H,oi (4.57) 

Where Hpot, Hsoi are the closure (with respect to the intrinsic norm \\.\\h) of the sets 
of Cpot, Csol the sets of smooth, T^-periodic, potential and solenoidal vector fields, i.e. 
with C = (C°^(T'^))"' 

Cpot = G C I 3/ G C°°(T'^) with C = V/} (4.58) 

Csol = e C\3p G C with div(p) = and e = pexp(2C/) j e-^^^^^dx} (4.59) 
Thus H is a real Hilbert space equipped with the scalar product 

{^■,^)h = j^'^ i{x).i'{x)mu{dx) and by the variational formulation (|2.2|), for I G W^, 
^ID(U)1 is the norm in H of the orthogonal projection of I on Hsoi and I = Vxi + 
exp(2C/)p/ is the orthogonal decomposition of /. 

^/HD(U)l = dist(/, Hpot) (4.60) 

Now by duality for all G -fT 

d\si{i,Hpot) = sup (4.61) 

From which we deduce the following variational formula for the effective diffusivity by 
choosing ^ = / G M*^ 

HD{U)1= sup Uj^^-Pdxf 

pec div(p)=o JjdP^ exp{2U)dx J^^ exp{-2U)dx 



Note that the equation ( 4.62| ) gives back Voigt-Reiss's inequality by choosing p = I. 



Let Q{U) be the positive, definite, symmetric matrix given by the variational formula 
( 4.55| ) . Then the following proposition is a direct consequence of the equation ( [4.62| ) . 



Proposition 4.3. For all I G S'^"^ 

ID{U)1 J^^Q^p(^2U)dx Jjaexp{-2U)dx ^lgF-i^(,Q{U)^ (4.63) 

Choosing an orthonormal basis diagonalizing Q(U), it is an easy exercise to use this 
proposition in order to establish a one to one correspondence between the eigenvalues of 



Q(U) and D(U) to obtain the proposition 4.2 



4.1.4.1 Dimension two In dimension two, the Poincare duality establishes a sim- 
ple correspondence between Q{U) and D{—U). 

Proposition 4.4. For d = 2, one has 

Q{U) = ^PD{-U)P (4.64) 
where P stands for the rotation matrix 

'1 ~o} (*"^> 
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Indeed by the Poincare duality one has 



sol 



{PVf : f G C°°(T"')} and the 



proposition 4.4 fohows from the definition of Q{U). The proposition [4.1| is then a direct 
consequence of the proposition AA and one deduces from the equation ( [4.53 ) that if 
D{U) = D{-U) then 



XmUDiU))Xrain{DiU)) 



ex.p{2U)dx / ex.p{-2U)dx) Id 



Which leads to the equation ( 3.17 ) of theorem by theorem 3.1 of |39|. 



(4.66) 



4.2 Sub-difFusive behavior from homogenization on infinitely 
many scales 

4.2.1 Anomalous behavior of the exit times: Theorems 13. 4] and 13.5. 



4.2.1.1 In this subsection we will prove the asymptotic anomalous behavior of the 
mean exit times E2,-[T(0,r)] defined as weak solutions of Lyf = — 1 with Dirichlet con- 
ditions on dB{0,r). Here C/„ S C^(T'^), nevertheless we will assume first that those 
functions are smooth and prove quantitative anomalous estimates on Ei;[r(0, r)] depend- 
ing only on the values of D(Vq^), Kq and Ka- Then, using standard estimates on the 
Green functions associated to divergence form elliptic operators (see for instance fl^ ) 
it is easy to check that the exit times E3;[r(0, r)] are continuous with respect to a per- 
turbation of V in L°°(i?(0, r))-norm. Using the density of smooth functions on 5(0, r) 
in the set of bounded functions, we will then deduce that our estimates are valid for 
Un e C"(T'^). 

Thus we can see the exit times as those associated to the solution of ( 3.23| ) and take 
advantage of the Ito formula. 



The central lemma of the proof is lemma |3.l| , which will be proven in the paragraph 
4.2.1.2| . 

Writing mlj(dx) = e~^^^^^ dx{J^^^ g"^^*^^') dx)~^, we will prove in the paragraph 4.2.1.3 
that for P e C°^(5(0,r) 



/ E^+^[r(0,r)]m[,+p(dx) < e^O-B.o.jCP) f ^'^ [r(0, r)] m[,+p(dx) 

JB{0,r) JB(0,r) 



(0,r) 

> g-2 0scB(o,.)(-P) f E^[T(B(0,r))]m^+p(da 

JB(0,r) 



(4.67) 



We give here the outline of the proof (see |39| for d = 1). A perpetual homogenization 
process takes place over the infinite number of scales 0, . . . , n, . . . and the idea is still 
to distinguish, when one tries to estimate ( 3.29^ , the smaller scales which have already 
been homogenized (0,... ,nef called effective scales), the bigger scales which have not 
had a visible influence on the diffusion {n^ri, ■ ■ ■ ,oo called drift scales because they will 
be replaced by a constant drift in the proof) and some intermediate scales that manifest 
the particular geometric structure of their associated potentials in the behavior of the 
diffusion (uej + 1,... ,nrfri — 1 = "-e/ + iT-per Called perturbation scales because they 
will enter in the proof as a perturbation of the homogenization process over the smaller 
scales). 

We will now use ( p5D and ( pT]) to prove theorem p.4| . For that purpose, we will first 
fix the number of scales that one can consider as homogenized (we write ef for effective) 



riefir) = sup{n > : g^^+^^^^^^^+^^^^^ii^ < C'i/(8Crf)r2} < oo 



(4.68) 



23 



where and Cd are the constants appearing in the left term of ( ^.25 ), next we fix the 
number of scales that will enter in the computation as a perturbation of the homoge- 
nization process (we write per for perturbation) 



n. 



.per(r-) = inf{n > : Rn+i > r} - ne/(r) (4.69) 



For r > Cd,Xo,pmax5 ^e/(^') and riperir) are well defined. Let us choose U = Vq^^^''\ 
^ = (EH)' ^6 ^ill bound from above, OscB(o,r)(K^^(^r)+i) 



'^^'^V Vj.(r)+1 + II V/(r)+nper{r)+lll« 



In the lower bound of ( [4.67| ) when x € B{0,r/2) we will bound E^[r(0,r)] from below 
by E^[r(j;, r/2)] and in the upper bound when x € -6(0, r) we will bound it from above 
by IE^[r(x,2r)]. Then using ( p. 251 ) to control those exit times one obtains 

2 

E^\T{B{0,r))]m^'^''''\dx) < Cde'^^-"+«"--W^° — ^-^ ^ 

Bio,r) V'^iJ y y }- d X^..{D{V^'-^f(-))) ^^^^^ 

Theorem |3.4| follows directly from the last inequalities by using the estimates (^]^) on 
-^(^")) ( |2.11| ) on Rn and observing that 

nperir) < inf{m > : ^^+^^f('^+' > c,eMr)+2)i9d+i5)Ko/2^ (4,71) 
The proof of theorem 3.5 follows similar lines, the stability result ( [4.67| ) being replaced 



by the stability condition |3J. 



4.2.1.2 It is sufficient to prove the equation ( 3.25| ) for x = 0. 

Let for I € S'^^^, xi be the T^-periodic solution of the cell problem associated to Lu 
with xz(0) = 0. 

Write the T^-periodic solution of the ergodicity problem = \l — Vx^P — ^ID{U)1 
with (/)/(0) = 0. Write Fi{x) = l.x — Xi{x) aiid Tpi{x) = Ff{x) — (j)i{x), observe that since 
LuF^ = \l- VxiP it follows that 



Lu^pi = HD{U)1 (4.72) 
The following inequality will be used to show that X^iLi V'e^ behaves like \x^ 

d 

CM"" - C2{\\x.\\lo + U.Woo) < j;V'e,(x) < Csdxp + llx.IlL + U.Woo) (4.73) 

i=l 

Using theorem [4.1| (theorem 5.4, chapter 5 of |^^) to control Fi and ifji over one period 
(observing that LijFi = 0, Lijil)i = —1) and using xi = ^-^ ~ Fi and (pi = Ff — ipi one 
obtains easily that H^. ||oo < Cde^^'^~^^^^ '^'^^^^^ R'^ , combining this estimate with ( 4.16 ) one 
obtains 

\\X. \\l + moo< Cde^'"^''^ (4.74) 
Since V has been assumed to be smooth (in a first step), we can use Ito formula to obtain 

i^i{yt)= fv^i{ys)du:s + HD{U)lt (4.75) 
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Now (ej being an orthonormal basis of M'^) write Mt the local martingale 

d 



Mt = Yl V'e. (yt) - Trace {D{U)) t (4.76) 



i=l 

Define 



r'(0,r)=inf{t>0 : \Y^^,^{yt)\ = r} 



i=l 



According to the inequality ( 4.73 ) one has 



r'(0,Cir2 - Cadlx.IlL + ll'A.lloo)) < r(0,r) < r'(0,C3(r2 + \\x.\\l + ||</..||oo)) (4.77) 



Since M(^T-/(o,r) is uniformly integrable (easy to prove by using the inequalities ( 1.77D ) 
one obtains 



Thus, by using the inequality ( [4.74| ) and the Voigt-Reiss' inequality D{U) > e 2 0sc((/) 
one obtains the equation ( ^.25 ). 



4.2.1.3 The proof of the weak stability result ( 4.67] ) is based on the following obvious 



lemma that describes the monotony of Green functions as quadratic forms, i.e. 

Lemma 4.8. Let Q be a smooth hounded open subset o/ M^. Assume that M,Q are 
symmetric smooth coercive matrices on 0,. Assume M < XQ with A > 0, then for all 
f € C^{Q), writing Gq the Green functions of —VQV with Dirichlet condition on d^l 



Ggix, y)f{y)f{x) dxdy<\ Gm{x, y)f{y)f{x) dx dy (4.79) 

Jn 

Proof. Let / G C^iTt). Write tpM , tpQ the solutions of -VAIV^jm = f and -S/QVipq = 
f with Dirichlet conditions on Observe that ipM and ipq are the unique minimizers 
of lM{h, f) and /q(/i, /) with 

lM{h,f) = l [^VhMVhdx- [ h{x)f{x)dx (4.80) 

and luii'M, f) — ~\ '^M{x)f{x) dx. Observe that since M < XQ 

lM{hJ)<XlQ{hJ-) (4.81) 



and the minimum of the right member in the equation (4.81) is reached at 4^q/X. It 



follows that ipQ{x)f{x) < X ipAf{x)f{x), which proves the lemma. □ 

Then, the equation (}4.67| ) follows directly from this lemma by choosing Q = e~'^^^^^\ 
M = e-"^^ and observing that E^[T(0,r)] = 2 /g^g^^^ G^-2u i^{x,y)e-'^^^y^ dy. 
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4.2.2 Anomalous heat kernel tail: theorem B.6 



4.2.2.1 From the pointwise anomaly of the hitting times of theorem 



one can 



deduce the anomalous heat kernel tail by adapting a strategy used by M.T. Barlow and 
R. Bass for the Sierpinski Carpet. This strategy is described in details in the proof of 
theorem 3.11 of |^ and we will give only the main lines of its adaptation. 
We will estimate ¥x[T{x,r) < t] and use Pa;[|2/t| >''']< ¥x[T{x,r) < t] to obtain theorem 



3.6 



Using the notations introduced in theorem and M := {d, Kq, Ka,a, /x, Amax)) it will 
be shown in paragraph 4.2.2.2 that for r > C(M, pmax) one has 

Now we will use lemma 3.14 of |^ given below (this is also lemma 1.1 of |^). 

Lemma 4.9. Let ^1,^2, • • • ,Cn, V be non-negative r.v. such that V > Y17=i ^i- Suppose 
that for some p G (0, 1), a > and t > 

IP(e^ <tk(6,-- - <P + at 

Then 

ln¥(V <t)< 2( — )^ -nln- 
p p 

Let n > 1 and g = ^. Define the stopping times Si i > hj Sq = and 

Si+i=mi{t>Si : \yt-ys,\>9} 

Write = Si — Si-i for i > 1 Let J-t be the filtration of yt and let Gi = J^Si Then it 
follows from the equation ( [4.82| ) that for 

r/n > C{M, p^ax) (4.83) 



¥.[C.+i<t\{G)i]=Fy,^[T{ys,,g)<t] 
Since — ySi+J = 5 it follows that a.s. \x — ys„ \ < r. Thus 



Sn = ^ii< T{x,r) 



i=l 

And by lemma |4.9| with 

a = C2i(M)(-)2+-('^)(i+^) p = l- C72o(M)(-)2-Mt 
r r 

One obtains 



lnP^[r(a;,r) < t] < 2 



1-C2o(?)2-M^ ) 



nln 



1 - C20(^)2-M^ 



(4.84) 



Minimizing the right term in ( 4.84| ) over n under the constraint (|4.83D and the assump- 
tions (|3|), Pmin > Cq^m, One obtains theorem pM. 
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4.2.2.2 The equation ( |4.82 ) is an adaptation of lemma 3.16 of [^]. Observe that 

^x[r{x,r)] <t + E,[l{T{x,r)>t)Ey,[T{x,r)-t]] 

<t + F^[l{T{x,r) >t)] sup Ey[T{x,r)] 

y(iB{x,r) 



Using Vy G B{x,r), Fy a.s. T{x,r) < T(y,2r) it follows by theorem 3.5 for r > 

C(M, pmax) 

C33(M)r2+'^W{i-T) < E^[T{x,r)] <t + F^[T{x,r) > t]C734(M)r2+'^W(^+^) 
Which leads to (148^ ). 
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